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We evaluate the effect of screening by bound electron in ( Be, e) +p — > ( S B, e) +7 transition in 
the framework of the adiabatic representation of the three particle problem. Comparison with two 
approaches (united nucleus and static) is presented. We discuss possible applications of this effect 
both for Solar Neutrinos and low energy fusion experiments. 
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I. INTRODUCTION 

o\ : 

0^ . In the recent years an increasing interest has been devoted to the accurate estimation of electron screening effects 
both for the solar plasma fusion rates and low energy Earth experiments (see Ref. and references therein). 

As a rule, one considers electrostatic screening in the solar plasma. This approach being classical or quantum 
correctly reflects the major properties of a process only for high relative velocity of the colliding nuclei, when electron 
density is supposed to be unchanged during the collision. 

In the case, when relative velocity of the nuclei is too small in comparison with the electron one, the electron density 
is changed accordingly to any relative configuration of the nuclei. Thus the electrons have an impact on a kinetic 
_ , energy shift of the nuclei. Such a process is considered in the adiabatic approach, which comes from the well known 
' Born-Oppenheimer method. 

An accurate treatment in the adiabatic approach of electrons from the continuum spectrum requires an additional 
research, while screening effects by bound electrons could be easily considered in the framework of the adiabatic 
approach (see, for example Ref. ||). 

For not very dense stars like the Sun, the commonly accepted Debye-Hiickel approximation on the calculation of 
screening effect is not quite adequate. Recently this question was discussed by A. Dar and G. Shaviv ||. This fact 
^\ • is usually of no importance considering that the screening due to the plasma electrons is by itself rather small (see 
""^ ' Ref. On the other hand, the lowdying bound electrons on a nucleus do screen the electric charge of the colliding 

^3 nuclei much effectively. Not only the stellar plasma but Earth experiments have such a feature. 
1 1 ' Indeed, laboratory experiments arc performed with atomic or molecular targets, and the ionic beam interacts with 
the target nuclei while the nuclei are surrounded by a number of bound electrons. It leads to that the fusion cross 
section is increased. At low and moderate energies, the fusion cross section of "bare" charged nuclei colliding with 
the relative momenta p and reduced mass M (in electron mass unit) is expressed as Ref. [ |l0| ]: 

I *{E) = ^§l e -^ (1) 

where S(E) is the so-called astrophysical factor which incorporates all nuclear aspects of the process, E is the 
collision energy of nuclei and r\ — M Z1Z2/ a e p( the latter factor comes from VWrf(O)- the Coulomb wave function of 
the internuclear motion at the origin) and a e is the hydrogen Bohr radius. 
The screened cross section differs by a factor j(E) defined as: 

1(E) , £~ = ^% (2) 

where a sc and G\, a are fusion cross sections both of the screened and "bare" nuclei, and V , b(0) is the wave function of 
the internuclear motion at the origin, taking into account the bound electron around. 

In this Letter we present the first quantum mechanical calculation of screening effect by bound electron in 

( 7 Be,e)+p^ ( 8 B,e)+ 7 

nuclear fusion. This reaction is of no importance for the Sun luminosity but it is of great interest in Solar Neutrino 
Puzzle. An analysis of Solar Neutrino Experiment Data (see Ref. shows that there is no room for 7 Be neutrinos. 
At the same time 8 B neutrinos are presented. Since S B nuclei result from the reaction: 

7 Be + p — > & B + ll (3) 
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it looks impossible to explain the existence of 8 B neutrinos and the absence of Be neutrinos in the frame of the 
Standard Solar Model. 

As one of the main consequences of electron screening effect applied to all solar fusions could be decrease in the Sun 
core temperature. Increased values of fusion cross sections could lead to the Sun cooling with the same observable 
Solar luminosity. However, just cool Sun models do not solve the Solar Neutrino Problem (see Ref. 

We show that bound electron essentially enhances the fusion rate in comparison with the reaction rate. Physics 
of this phenomena could be easily understood in the framework of united nucleus approach (see below). "Exact" 
solution of our problem is obtained in adiabatic approach for three particle problem. We compare our solution with 
the two relevant approaches (united nucleus and static approaches, described below) which give upper and lower 
estimate for the screening effect, in the considered here three-particle picture. 

A finite 7 Be — p interaction radius also is of great importance. A larger radius results in higher tunneling probability 
into the internal region and thus to a higher 7 Be{p, j) 8 B cross section. Theoretical model Ref. [|l3| also predicts 
increased astrophysical Si7(0)-factor. All above mentioned effects result in the increased 8 B production rate and the 
latter could be a really alternative to 7 Be + e — > 7 Li + v e reaction. 

II. METHOD OF CALCULATION 

We treat the Coulomb problem for three particles in the framework of the adiabatic representation. The basic 
two-center eigenfunctions are derived from the Schrodinger equation for two nuclei with electric charges Z\ and Z 2 
(in electron charge unit) (Z\ > Z 2 ) fixed on a distance R and for an electron around them. The Schrodinger equation 
is transformed into a infinite system of equations with separated variables. Our approximation consists in that we 
use only one two-center eigenfunction corresponding to the ground state of the system. There are some arguments 
for this approach. 

At first, the high energy states corrections (at fixed R) are of order of magnitude about the ratio of the electron 
mass to the proton mass. 

Then, excited energy levels correspond to the less energy of the united nucleus. It leads to their exponentially small 
contribution into the nuclear fusion rate in comparison with the ground state of the system. 

At least, only the ground state energy of the electron in the field of 7 Be and p nuclei, called ISa therm (see, for 
example Ref. 114]]) has the correct asymptotic behaviour (see below). 

The three particle wave function is presented as: 

^(R,f) = <j>(R)-^(R,r) 1 

where </>(i?) is the wave function of two colliding nuclei and ip(R,r) is the electron wave function, which depends on 
the internuclear distance R. 

The electron energy eigenvalue U n im(R) also depends on the internuclear distance. We use the tabulated values of 
Unim from Ref. (l^] for our purposes. 

In Fig. |l] the values for the ground state are plotted. At R — » the electron energy approaches to the energy of 

the united ion: U — ► — ii±l±£lL (energy unit here is 27.21 eV ) and for R — > 00 the electron energy approaches to the 

z 2 

energy of the isolated ion eZ\. U — ► — 

We consider the case Z\ = 4 and Z 2 = 1- In the adiabatic approach the electron energy U n i m (R) in the field of 
two nuclei serves as an effective attraction potential for the nuclei. Evidently the main effect of the electron screening 
comes from the collisions with zero orbital moment of nuclei. The latter and the correct behaviour of U n i m (R) (at 
R — > and R — ► 00) for the ground state allow us to consider the case n = I = m = 0. Let us denote t/ooo as U(R). 
With the electron energy U(R) in hand one can solve the scattering problem for two nuclei 7 Be and p. The total 
effective potential reads: 

V(R) = ^ + U(R) (4) 

We calculate values of the wave function of relative motion of the two nuclei at kinetic energies 0.1 keV < E < 100 
keV and compare the calculated ^(O)! 2 with |^ , co «;(0)| 2 = e ^Z v _ 1 - value of the scattering wave function of two 
particles in the origin. 

For simple estimations one can use well known united nucleus approach that consists in following. Consider a 
fusion of two nuclei from initial states with bound electrons. Let the total negative energy of these bound electrons 
be E1+E2. In the final state there is an united nucleus with bounded electrons around and let the total energy of the 
electrons in this new ion be E u . Thus the difference AE = E\ + E 2 — E u adds to the kinetic energy of the two nuclei 
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accelerating them. So, the united nucleus approach changes | 'I' £■ (0) | 2 to e+ae(0)\ 2 ■ This is a rather good approach 
since the main contribution to the value |4 , £;(0)| 2 comes from internuclear distances less then classical return point 
in the potential. 

Since united nucleus approach replaces decreasing E(R) (see in Fig. |l|) by the constant AE this approach serves as 
an upper estimate. 

Another approach uses static wave function of the bound electron. It means that this wave function does not 
depend on the internuclear distance. Total three particle wave function is presented as a product: 

where the wave function for the electron bounded on the nucleus with electric charge Z\ reads: i\) s {f) = y ir , 

and 4> S (R) is the wave function of two nuclei. Averaging the electron coordinates , 4> S {R) can be derived from the 
Schrodinger equation with the potential: 

V st (R) = ^-2k- (^~~~J[~~ ~ Z ^ 2RZl ) (6) 

Kinetic energy of the colliding nuclei is counted from the electron eigenvalue energy in the field of Z\ nucleus. This 
approach works well for nuclei with relative velocity v much greater then electron velocity v e : v >> v e . 

Attractive part of the potential V s t approaches to the value \Z\Z-i (in 27.21 eV unit) at R — ► 0. Since Z\Zi < 

^ Zl+ 2 Z2 ^ — the value of the electron energy in united ion, the static approach gives a lower estimate for an enhancement 
factor. 

We use all three approximations in order to compare their validity. 

III. RESULTS 

The numerical solution of the scattering problem of the two nuclei in the potential (Q) was obtained on the mesh on 
R for R = and R = R ma x in Numerov's scheme. We varied step h and R m ax in order to ensure that final result is 
not changed substantially. Also for checking purposes we reproduced the value of \^ CO ui{^)\ 2 substituting U{R) = 0. 

The enhancement factor (^) is plotted in Fig. g. In Fig. ^ we compare the "exact" numerical solution with both 
united nucleus and static approaches. 

At kinetic energies above 2 keV a very good agreement between all three approaches is obtained, though at lower 
energies the united nucleus approach overestimates and the static approach underestimates an electron screening 
effect. It is easy to see that simple united nucleus prescription gives much closer values for electron screening effect 
then the static approach. Thus the united nucleus approach being compared with "exact" numerical solution could 
be easily used in stellar evolution code. 

As one can see, electron screening could essentially increase nuclear fusion rate and must be taken into account 
both for Earth experiments at low energies and for an understanding of stellar phenomena. 

Deficit of 7 Be neutrinos in comparison with 8 B neutrinos (called "second Solar Neutrino Problem", see Ref. p6|) 
could be a good example of this effect. At the same Solar Model input data the flux of 8 B neutrinos increases due to 
the increased rate of 8 B nuclei production and the flux of 7 Be neutrinos decreases due to the alternative reaction (^). 
At mean kinetic energy in the Sun core at temperature T = 15Tq (Tg = 10 6 K) this effect is 15 times of magnitude. 

However due to the Gamov's "window" around 17 — 20 keV the screening effect of the bound electron is reduced to 
10%. But if some nuclear fusion reaction go faster in the Sun core (due to the electron screening) then the temperature 
in the Sun core could be lower. Under this assumption it might be possible to explain the Solar Neutrino Problem. 

In conclusion, we have performed quantum mechanical calculation of the screening effect taking place in ( 7 Be, e) + 
p — > ( 8 -B, e) + 7 reaction at finite kinetic energy of the nuclei. We compared our results with two approaches giving 
upper and lower limits for the screening effect and studied the validity of the approaches. Possible application to the 
Solar Neutrinos is briefly discussed. 
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FIG. 3. j(E) factor for united nucleus (dotted line), "exact" numerical solution (dashed line) and static approaches (solid 
line) 
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